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Nuclear matter equations of state based on Skyrme, Myers-Swiatecki and Tondeur interactions 
are written as polynomials of the cubic root of density, with coefhcients that are functions of the 
relative neutron excess 5. In the extrapolation toward states far away from the standard one, it 
is shown that the asymmetry dependence of the critical point (pc,5c) depends on the model used. 
Ifowever, when the equations of state are fitted to the same standard state, the value of 5c is almost 
the same in Skyrme and in Myers-Swiatecki interactions, while is much lower in Tondeur interaction. 
Furthermore, Sc. does not depend sensitively on the choice of the parameter 7 in Skyrme interaction. 

PACS numbers: 21.65.+f, 21.30.Fe 



I. INTRODUCTION 

Nuclear matter is considered as an uncharged nucleon system distributed uniformly in the space, and nuclear matter 
equation of state is the energy per nucleon e(/3, 5) of nuclear matter given as function of nucleon density p and relative 
neutron excess 5. The equation of state e(/3, 5) is a fundamental quantity in theories of neutron stars and supernova 
explosions, as well as in studies of nucleus-nucleus collisions at energies where nuclear compressibility comes into play 

The main measured quantities which can provide information about equation of state (EOS) are the binding energies 
and other data of finite nuclei. As the finite nuclei are in states near the standard nuclear matter state, which is 
the ground state of nuclear matter with normal nucleon density po and zero neutron excess, 5 — Q, therefore, our 
knowledge about EOS can be confirmed experimentally only in a small region around p ^ and (5 0. In this 
region, the main quantities which specify the EOS are the coefficients oi (volume energy), J (symmetry energy), Kq 
(incompressibility), L (density symmetry), and Kg (symmetry incompressibility) . Nowadays the quantities which are 
known with enough precision are oi, J and Kq, while the last two are still under investigation. 

However, there is currently considerable interest in the very neutron rich nuclei and the energetic heavy-ion collisions 
where the nuclear matter state is beyond this region. As any direct information beyond this region is difficult to come 
by, extrapolation is inescapable and, in this case, a nuclear model is required. This model is fitted to binding energies 
and other data of finite nuclei at first, then applied to nuclear matter to derive the EOS. In this way, the obtained 
EOS can be considered as being fitted indirectly to a region around the standard state, but its prediction on states 
beyond this region should be regarded as an extrapolation. Obviously, the reliability of this extrapolation depends 
on the foundation of the model. In order to be reliable, the model should be based on a well-founded theory with 
as few adjustable parameters as possible, which are fitted to as many high accuracy measured data as possible. 
Basically, this is what we understand by model of effective nucleon-nucleon interaction, for example the Skyrme , 
the Myers-Swiatecki and the Tondeur interactions jj]. The EOS's derived from these effective interactions have 
analytical expressions which can be pictured and calculated easily, so they are used widely in the literature, even 
they are based on energy functional theories which are in the macroscopic level. Instead of effective nucleon-nucleon 
interaction, nuclear matter is studied also in more basic level by "microscopic" potentials available in the market and 
by sophisticated many-body theories for many years Q Q ||] f^- 

The purpose of the present paper is to discuss the EOS's given by Skyrme ||], Myers-Swiatecki ||] and Tondeur 
interactions Q, in comparison with the microscopic calculations. In Section 2, the nuclear EOS given by these 
interactions is presented. The equilibrium condition and the properties of standard nuclear matter are discussed in 
Sec. 3, while the predictions for nuclear matter away from the standard state are given in Sec. 4. Sec. 5 makes 
comparison with some microscopic calculations. In Sec. 6 a short discussion and summary are addressed. Appendix 
A presents a specific discussion on Myers-Swiatecki interaction and Appendix B gives some formulas used in Sec. 4 
to calculate the interaction parameters from standard nuclear matter quantities. 
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II. NUCLEAR EQUATION OF STATE 



The nuclear energy Ejsi of a nucleus can be written as 

En = j d\£N, (1) 

where the nuclear energy density functional £n can be written with enough generality as 

£n ^ p{v)e{p,5)+£GD. (2) 

where 

£gd - \Qi{^pf + Q2 [(Vp„)2 + (Vpp)2] (3) 

is the gradient density dependent term. In the above equation, p ^ Pn + Pp, Pn and pp are the neutron and proton 
densities respectively, 

S = ^^^^ (4) 
P 

the relative neutron excess, Qi and Q2 the model parameters related to the finite size and the surface effects of nuclei. 

The EOS e(n,S) depends on the model of interaction, while the functional Sgd depends also on the model of 
nuclei. Eq. is exact for Skyrme and Tondeur interactions, and is approximate for Myers- Swiatecki interaction 
(see Appendix A). However, the functional Sgd is irrelevant to the present discussion, since it is irrelevant to the 
nuclear matter property. In the following the EOS's based on Skyrme, Myers- Swiatecki and Tondeur interactions will 
be given. 



A. Skyrme interaction 



The EOS based on Skyrme interaction can be written as 
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where po = 3/47rrQ, vq the nuclear radius constant, T an appropriate constant with dimension of energy such that 
the D coefficients are dimensionless, and 7 a model parameter. It is convenient to choose T as the Fermi energy of 
standard nuclear matter. 
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where m is the nucleon mass. The D coefficients are 
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and to, ti, t2, t3, xq, xi, X2, X3, 7 are the interaction parameters. It is worthwhile to note that, among these interaction 
parameters, only to, t^, xq, X3, Si, S2, and 7 appear in the EOS and thus are relevant to the nuclear matter properties. 
Beside these interaction parameters, there is another interaction parameter Wq ^ that appears only in the coefficients 
Qi and Q2 and thus is irrelevant to the EOS. 
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B. Myers-Swiatecki interaction 



The EOS based on Myers-Swiatecki interaction can be written as 



where 
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Dl'HS) = ^{Bi [(1 + 5f/^ + (1 - 5)«/3] + i?„(l - S') [(1 + <5)2/3 + (1 _ 5)2/3] I 



(15) 



In the above equations, a, Bi, B^, "fi, 7u, and ^ are the interaction parameters. In addition to these parameters, there 
is another one a the Yukawa range of force, that is irrelevant to the EOS, as it appears only in the coefficients Qi 
and Q2. 



C. Tondeur interaction 



The EOS based on Tondeur interaction can be written as 
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In the above equations, a, 6, c, and 7 are the interaction parameters. In addition, there are another two interaction 
parameters d and rj M, that are irrelevant to the present discussion as they appear only in Qi and Q2- 



III. STANDARD NUCLEAR MATTER 

The EOS given in the last section can be written generally as 
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The equilibrium condition de/dp\Q ~ 0, by which the standard state p = pg at 5 = is defined, gives the following 
relationship among 152(0), £'3(0), -05(0), £'7(0) and 7: 



202(0) - SDsiO) + 5D5(0) + 7^7(0) = 0. 



(21) 
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The 5 quantities of nuclear matter ai, Kq, J, L, and can be expressed as: 



ai = -e(po, 0) = -| [i?20 - 27^50 - (7 " 3)D^o] , 



(22) 
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= T[- 2D22 + IOD52 + 7(7 - i)0^2] , 



where 



Relation ( |2lD is used in obtaining Eqs. (^2|) and (|23|), from which the following formulas can be derived: 

Ko = 15ai + [3O20 + (7 - 5) (7 - 3)i?^o]r, 



Ao = A(0), A2 = ^^^|o, i = 2,3,5,7. 



(26) 
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(28) 



Ko = 3701 + [(7 - 2)020 - 2(7 - 5) Ao]T. (29) 
The specific discussion for Skyrme, Myers-Swiatecki and Tondeur interactions will be given in the following. 

A. Skyrme interaction 

For Skyrme interaction, the following relationship can be obtained from Eq. (pT|): 



8T16 T \ 2 J T ^ 2 ' 5 



(30) 



Therefore, among 7 parameters to, t^, xq, 0:3, si, S2, and 7, only 6 of them are free. Considering this relation, it can 
be shown that oi, Kq, J, L and Ks are independent each other, in the Skyrme EOS. 
A relation connecting to ai and Kq can be obtained from Eq. (^8|), 
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If t3 = 0, as fli \QMeV and T ~ 37A/ey are well-known from measurements, this formula gives the estimation 
Ko ^ 306A/ey. Therefore, in order to have Kq value lower than 306Mel^, the fourth term (p/po)'*^^ in the Skyrme 
EOS is needed. 



B. Myers-Swiatecki interaction 

For Myers-Swiatecki interaction, 0^-'^{6) = 0, all the 7-dependent terms in Eqs. (|2l])-(|2^) do not appear. The 
equilibrium condition (21) can be transformed into the following relation among a, B and 7: 



5a - 105-4(1-7) = 0, 
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where B and 7 are defined respectively as 

Bi,u^l{lT(:)B, 7/,« = ^(1tC)7- (33) 

Therefore, there are only 4 independent interaction parameters in the Myers-Swiatecki EOS, a, B, ^ and C, if 7 is 
solved from Eq. ( |3^ ) as a function of a and B. Correspondingly, there are only 4 independent variables among ai, 
Kq, J, L and Kg in the Myers-Swiatecki EOS. Actually, the following relationship can be derived: 
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Furthermore, formula ( p8| ) for the Myers-Swiatecki EOS becomes 

A-o = 15ai + p(l-7)T. (37) 
5 

For 7 = 0, Myers-Swiatecki interaction is reduced to Seyler-Blanchard interaction and this formula gives the 
estimation Kq ~ SOQMeV, the same as that discussed for Skyrme interaction. Hence, 7-dependent terms in Myers- 
Swiatecki EOS are required, in order to obtain Kq lower than 306MeV ]lO| . 

C. Tondeur interaction 



For Tondeur interaction, D^{S) — 0, the term involving I?5(0) in Eq. (21) as well as all the terms involving D 



^50 



and Z?52 in Eqs. (|22|)-(|26|) do not appear. The equilibrium condition (^IJ) now is a relation among a, b and 7: 

^ + 2^ + 6=0. (38) 

In this case, there are only 3 independent interaction parameters in Tondeur EOS, for example a, c and 7. Corre- 
spondingly, there are only 3 free variables in ai, Kq, J, L and Kg, for example ai, Kq and J, since it can be shown 
that 

L = 2J, = -2J. (39) 
In addition, the following relationship among Kq, ai and 7 can be written for Tondeur EOS from Eq. (|29|): 

= 3701 + 1(7 -2)r. (40) 
5 

From ai ~ IGMeV.T ^ 37MeV and i^o ^ 220MeV, it can be evaluated that the appropriate integer is 7 = 4, as 
given by Tondeur |4|. In this case, i.e., if 7 = 4 is chosen, there are only two interaction parameters to be freely 
adjusted in the data fit, for example a and c. Correspondingly, there are only two independent variables in ai, Kq, 
J, L and Kg, for example oi and J, when Kq is calculated by Eq. (^). From oi ~ IGMeV, T ~ 2,7 MeV and 7 = 4 
we can evaluate Kq ~ 236MeV. It is worthwhile to note that the value given by Tondeur is Kq = 235.8Mey 

The equilibrium condition is checked by calculating the expression on the lefthand side of Eq. ( ^l|) for Skyrme, 
Myers-Swiatecki and Tondeur interactions, using the interaction parameters and the nuclear radius constant tq given 
in Refs. and Q, respectively. These parameters will be referred to as the original interaction parameters 



thereafter. Besides, the following physical constants |13| are used in the present calculation: he = 197.32891MeT^- /m, 
= 938.90595Mey/c^ 
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The calculated values are given as EC in the second column of Table 1. It shows that the equilibrium condition 
of standard nuclear matter is fulfilled in the data fit to determine the original parameters of Skyrme (lst-5th row), 
Myers-Swiatecki (6th row), and Tondeur interactions (7th row), respectively. 

The standard nuclear matter properties ai, ifg, J, L, and Ks, calculated from Skyrme (lst-5th row), Myers- 
Swiatecki (6th row), and Tondeur interactions (7th row) respectively, are also given in the 5th-9th column of Table 
1, all in MeV. In this table, 7 is a model parameter in Eqs. (||) and ([l6| ) for Skyrme and Tondeur interactions 
respectively, ro the nuclear radius constant used in the respective interaction, in /to. 

As a comparison, the last two rows of Table 1 (labeled by CWS) present the result obtained by fitting ai, Kq, J, 
L, and Kg directly to nuclear masses It can be seen that these quantities have values close each other, except 
the case SIII, where the value of Kq and Ks is far away from others. The average over the 2nd to 7th row gives 
ai = 15.97Mey, Kq = 234.4Mey, J = 29.25Afey, L = 48.63Mey, and = -126.9MeV. 



IV. NUCLEAR MATTER AWAY FROM THE STANDARD STATE 



The nuclear matter state with zero pressure and minimum energy per nucleon can be solved from the following 
equation: 



de 
dp 



= 0. 



(41) 



Usually there are several solutions, we should choose that one has minimum energy per nucleon. This solution gives 
density as function of 6: 



Pm = Pm{S). 

For (5 = 0, Eq. ( ^ ) is reduced to the equilibrium condition of standard nuclear matter, we have 
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The incompressibility of non-equilibrium nuclear matter, which is of interest in many applications, can be defined 
as M 



Kip, 6)^9 
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where P — p^de/dp is the pressure. Along the line of minimum (p^), this K{p, 5) becomes 



(44) 
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At the standard state (po;0) we have K^ijS) = Kq. At the critical point {pc,^c)i where the maximum and the 
minimum are coincident, the curvature of e{p,5c) versus p changes sign and Kmi^c) = 0. So Km{5) starts with 
K{j and ends at when 5 increases along the line of minimum. In addition, the generalized symmetry energy of 
non-equilibrium nuclear matter can be defined as § Q |] 
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In term of this quantity, the usual symmetry energy J can be expressed as 

J = Jipo)- 



(46) 
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For nuclear matter not far away from the standard state (po,0), the EOS can be written approximately as |14 

1 „oN/iO-/9n\2 Lfp-po 



In this approximation, we have 
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Using Eg. (^8|) , the following solutions can be obtained: 



J(p)-J+^^+^f^)^ (50) 



Pm^Po{l-^f), (51) 

e„ « -ai + JS"^, (52) 

Km-Ko + KsS^, (53) 

where only the linear term in S"^ is kept. The systematics of nuclear central densities based on elastic electron 
scattering data [^t) jl^ and muonic atom spectroscopy data [|l8| provide a direct evidence for Eq. (|5l|). 

Thus, in the plot e{p,d) versus p, we have the geometric meaning of oi, Kq, J, L, and Kg'- the standard state is 
at the minimum point pm = po with depth ai and curvature proportional to Kq; when the minimum is moved with 
increasing S from 0, the decrease of pm is controlled by 3L/Ko, the increase of depth is controlled by J, while the 
decrease of curvature is controlled by —Kg. Therefore, the quantities ai, Kq, J, L, and Ks are characteristics of 
nuclear matter not only at standard state but also at the state not far away from the standard one. In this way, the 
interaction with different value of these quantities will predict different properties of nuclear matter that are not far 
away from the standard state. 

The exact solution Pm{5) depends on the interaction. The analytic solution is possible for SIII, Ska, Myers- Swiatecki, 
and Tondeur interactions, while the numerical solution is appropriate for SkM, SkM*, and RAPTP interactions. 

For Myers-Swiatecki interaction, (^ij) is a cubic equation which gives 



'""^ =2sos*«(g + 3), (54) 
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The superscript AIS for Myers-Swiatecki's D is dropped for simplicity. The critical 5c, where the maximum and the 
minimum is coincident, is determined by 



D2{&c) = 5sgZ?5('^c), (56) 



which corresponds to 6 = Q and 



Pc = Pni{5c) ^ sIpq. (57) 
For Tondeur interaction with 7 = 4, ( pl| ) is a quadratic equation which gives 

= ^i^^^ + [9^3 - i2DM6)]'/'}, (58) 

\ Po y 0-L>i 

where D3 and D4 are numbers, the superscript To for Tondcur's D is dropped also. The critical point is given by 

9Dl - S2DMSc) = 0, pc= (^)%o - -{^Y- (59) 

As the location po and depth oi are different for different equation of state, as shown in Table 1, a way to make 
comparison is to plot the normalized energy per nucleon e/oi as a function of the relative nucleon density p/po for 
given 6. Fig. la shows this e/ai versus p/po for 5 = 0, calculated by the various Skyrme interactions (solid lines), 
Myers-Swiatecki interaction (dot-dashed line), and Tondeur interaction (dashed line). The solid lines in the righthand 
side of the plot, from top to bottom, correspond to SIII, Ska, RATP, SkM, and SkM* interactions respectively. The 
difference between SkM and SkM* is negligible and Myers-Swiatecki is almost coincident with RATP. This sequence 
is just the decreasing sequence of i^o's value, as shown in Table 1: the smaller value of Ka, the smaller curvature of 
the curve at the standard state, thus the softer the EOS. 
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A natural question is: what is the difference among these EOS's, if the standard state is the same with same 
location po, depth ai, curvature ~ Kq and so on? In order to make this comparison, the interaction parameters 
should be readjusted according to chosen po, ai, Kq and so on. For the value of po, we choose tq = 1.140/to which 
is well determined by the data fit to nuclear charge radii jl^] extracted from elastic electron scattering data [0. In 
addition, we can choose the value of ai, Kq, J, L, and Kg m an appropriate way. In this case, Eq. (j^) should be 
fulfilled for Myers- Swiatecki interaction, while Eqs. (|3^) and (^ ) should be fulfilled for Tondcur interaction. The 
chosen values used to calculate the interaction parameters are listed in Table 2, while the formulas used to perform 
this calculation are given in Appendix B. Among these values, ai and J are the average values given in the last 
section, Kq and L are calculated by Eqs. ( |40| ) and (39). Kg is calculated by Eq. (|3|) in Tondeur's case while by Eq. 
( ^ ) in Myers-Swiatecki's case. In Skyrme's case. Kg can be chosen from either Myers-Swiatecki's or Tondeur's value, 
there is no significant difference in the calculated result which will be shown in the following. 

The calculated Skyrme interaction parameters are given in Table 3. The Myers-Swiatecki interaction parameters 
are calculated as: 



a = 2.06285, 5 = 1.05232, 7 = 1.05222, ^ = 0.12333, C = 0.37363. (60) 

For 7 = 4, Tondeur interaction parameters are calculated as: 

a = -672.13Mey/TO^ b = 799.71MeVfm* , c = 99.n6MeVfm^ . (61) 

These parameters will be referred to as the readjusted interaction parameters thereafter. 

As a comparison with the result calculated by original interaction parameters, Fig. lb plots e versus p/ po for 5 = 0, 
calculated by readjusted interaction parameters. It can be seen that now there is almost no difference among these 
EOS's for 0.4 < p/po < 1.6. 

Fig. 2 displays in (a) the normalized nuclear incompressibility K/Kq versus the relative nucleon density p/po 
calculated by original interaction parameters at (5 = 0, and in (b) K versus p/ po calculated by readjusted interaction 
parameters at (5 = of the various Skyrme interactions (solid lines), Myers-Swiatecki interaction (dot-dashed line), 
and Tondeur interaction (dashed line). Considering the righthand side of the plot, in (a), the solid lines, from top to 
bottom, correspond to SIII, Ska, RATP, SkM, and SkM* interactions respectively. The difference between SkM and 
SkM* is negligible. In (b), the sohd lines, from top to bottom, are due to SkM, SkM*, RATP, Ska, and SIII, where 
SkM* is identical to SkM. Tondeur is coincident with Ska; SkM and RATP almost overlap. It can be seen from this 
figure that the difference among these curves is negligible for p/po < 1.2. This is expected from Eq. ( ^9[ ) which shows 
that the curve is determined essentially by Kq. 

Fig. 3 depicts in (a) the normalized symmetry energy J{p)/ J versus the relative nucleon density p/ po calculated by 
original interaction parameters, in (b) J(p) versus p/ po calculated by readjusted interaction parameters of the various 
Skyrme interactions (solid lines), Myers-Swiatecki interaction (dot-dashed line), and Tondeur interaction (dashed 
line). Considering the righthand side of the plot, in (a), the solid lines, from top to bottom, correspond to Ska, SkM, 
SkM*, RATP, and SIII interactions respectively. The difference between SkM and SkM* for p/po < 1.3 is negligible. 
In addition, the Myers-swiatecki's is almost coincident with that of SkM's. In (b), the difference among these curves 
is negligible. This is expected from Eq. (pO|) , which shows that, for the density p is not far away from po, the symmetry 
energy J{p) is determined essentially by J, L and Kg, and these quantities (J and L) are the same or almost the 
same {Kg) for the readjusted interaction parameters. 

Using the interaction parameters, we can calculate pm{5), em{5), and Km{5) along the equilibrium line. The result 
is shown in Figs. 4-6, while the critical point value (pc, 5^) is listed in Table 4, for Skyrme, Myers-Swiatecki, and 
Tondeur interactions respectively. 

Fig. 4 shows the equilibrium density pm as function of the relative neutron excess (5, calculated by (a) original 
parameters and (b) readjusted parameters of the various Skyrme interactions (solid lines), Myers-Swiatecki interaction 
(dot-dashed line), and Tondeur interaction (dashed line). The solid lines, in the middle range of 5, from top to bottom, 
in (a), correspond to SIII, RATP, SkM, SkM*, and Ska interactions respectively, the difference between SkM and 
SkM* is very small. In (b), the solid lines from top to bottom correspond to Ska, RATP, SkM, SkM*, and SIII 
interactions respectively. SkM and SkM* are the same whereas RATP, SkM and Myers-Swiatecki almost overlap. 

Fig. 5 gives the equilibrium energy per nucleon as function of the relative neutron excess 5, calculated by (a) 
original parameters and (b) readjusted parameters of the various Skyrme interactions (solid lines), Myers-Swiatecki 
interaction (dot-dashed line), and Tondeur interaction (dashed line). The solid lines, in the righthand side, in (a), 
correspond to RATP, SkM, SkM*, SIII, and Ska interactions, from top to bottom, respectively. SIII, Ska and Tondeur 
almost overlap, whereas the difference between SkM and SkM* is negligible. In (b), the solid lines from top to bottom 
correspond to SIII, SkM, SkM*, RATP, and Ska interactions, respectively, where SkM and SkM* are the same. SkM, 
RATP and Ska are almost coincident. 
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Fig. 6 plots the equilibrium incompressibility Km as function of the relative neutron excess 5, calculated by (a) 
original parameters and (b) readjusted parameters of the various Skyrme interactions (solid lines), Myers-Swiatecki 
interaction (dot-dashed line), and Tondeur interaction (dashed line). The solid lines in (a) correspond to SIII, RATP, 
Ska, SkM*, and SkM interactions, from top to bottom in the middle range of S, respectively. Ska, SkM* and SkM 
almost overlap. In (b), the solid lines, from top to bottom, correspond to Ska, RATP, SkM, SkM*, and SIII interactions 
respectively; SkM and SkM* are the same. 

V. COMPARISON WITH MICROSCOPIC CALCULATIONS 

In order to provide additional elements about the confidence on the effective interactions discussed above, it is 
interesting to make a comparison with some microscopic calculations which are based on a more fundamental level of 
theories as well as on very different physical input. In Fig. 7 the present predictions for the pure neutron matter EOS 
are compared with the theoretical estimates of Friedman and Pandharipande |^ , obtained from a variational framework 
based on the Urbana vu two-nucleon potential plus three-nucleon interaction model of Lagaris and Pandharipande 
pof . The neutron matter EOS e(p, 1) versus nucleon density p is calculated by (a) original parameters and (b) 
readjusted parameters of the various Skyrme interactions (solid lines), Myers-Swiatecki interaction (dot-dashed line), 
and Tondeur interaction (dashed line). Black diamonds denote the data taken from Ref. Q. The solid lines, in the 
right hand side of (a), correspond to Ska, SkM, SkM*, RATP, and SIII from top to bottom. In (b), the solid lines, 
from top to bottom, correspond to SIII, SkM, SkM*, RATP, and Ska, where SkM and SkM* are identical. It seems 
that in the low density region the neutron matter EOS's calculated by original parameters are closer to microscopic 
results than those by readjusted parameters. However, the situation is different if the density is extended to higher 
region (see Fig. 9). 

It is worthwhile to note that the generalized Skyrme interaction FPS21 proposed by Pethick, Ravenhall and Lorenz 
[0 has the property that it is a good fit to both the nuclear and neutron matter calculations of Friedman and 
Pandharipande. In this sense. Fig. 7 may be regarded also as a comparison beetween our results and those of FPS21. 

Even the EOS's based on effective interactions and energy functional theories discussed in the present work are 
essentially nonrelativistic, it is still interesting to see how they behave in the high density region. Fig. 8 gives the 
symmetric nuclear matter EOS's e{p,0) up to about lOpo, calculated by (a) original parameters and (b) readjusted 
parameters of the various Skyrme interactions (solid lines) , Myers-Swiatecki interaction (dot-dashed line) , and Tondeur 
interaction (dashed line). The full dots stand for the results taken from Ref. which is a microscopic calculation 
of EOS for dense nuclear and neutron matter based on the Argonne vn two-nucleon potential plus Urbana VII 
three-nucleon potential. The crosses denote the results taken from Ref. ||], which studied the properties of dense 
nucleon matter and the structure of neutron stars, using variational chain summation methods and the new Argonne 
Vis two-nucleon interaction and the Urbana model IX of three-nucleon interaction as well as the relativistic boost 
correction to the two-nucleon interaction. The solid lines, in the right hand side of (a), correspond to SIII, Ska, RATP, 
SkM, and SkM*, from top to bottom, where SkM and SkM* are identical. In (b), the solid lines, from top to bottom, 
correspond to SkM, SkM*, RATP, Ska, and SIII, where SkM and SkM* are identical. The Tondeur's is very close to 
the Ska's. 

Fig. 9 is the same as Fig. 8 but for pure neutron matter EOS e{p, 1). The solid lines, in the right hand side of (a), 
correspond to Ska, RATP, SkM, SkM*, and SIII, from top to bottom. In (b), the solid lines, from top to bottom, 
correspond to SIII, SkM, SkM*, RATP, and Ska, where SkM and SkM* are identical. 

VI. DISCUSSION AND SUMMARY 

A discussion of the nuclear matter EOS's based on Skyrme, Myers-Swiatecki and Tondeur interactions is given 
in this paper. The equations are in the form of polynomials in the cubic root of density, with coefficients that are 
functions of the relative neutron excess and depend on the model of interaction. 

Most of the discussion about the nuclear EOS, up to now, focus at states around standard state, i.e. about the 
quantities ai, J, L, Kq, and Kg] especially Kq in supernova explosion and neutron star calculations and in heavy 
ion collisions. However, even these quantities or equivalently the interaction parameters were well-determined by the 
measured data of nuclei, mainly the nuclear masses, the extrapolation to states far away from standard state is still an 
open problem. It is seen that the difference among these EOS's is not significant in most of the relative neutron excess 
range which is of interest for both heavy ion collisions and supernova explosion calculations. However, if the equations 
are fitted to the same standard state, the equation based on Tondeur interaction is softer than others provided the 
relative neutron excess is not close to Pql . 



9 



The numerical result given in Section IV shows that the asymmetry dependence of the critical point depends on 
the model used in the extrapolation. When the EOS is fitted to same standard state, the Skyrme's and the Myers- 
Swiatecki's 6c are close each other, especially 5c does not depend sensitively on the choice of 7 in Skyrme interaction. 
On the other hand, the Tondeur's 5c is much smaller than others. This is because the value of the Tondeur's 5c 
depends sensitively on the interaction parameters, as it can be seen and checked numerically from the first equation 
of (I59I). In this content, in order to make a choice among these interactions for the extrapolation, experiments which 
can provide direct or even indirect information about nuclear matter with large asymmetry 5 and low density p are 
required. 
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APPENDIX A: 



It will be shown here that the energy density functional £gd of Myers- Swiatccki interaction can be written approx- 
imately in the form of Eq. (3h. In the Thomas- Fermi model of nuclei and up to the second order of the localized 
approximation given in Ref . [ I0[ , we have 



£'JE = ah{r/a)F('\r) + ^I,(r/a)F('Hr), 



(Al) 



where a is the Yukawa range of force. 



X 



(A2) 



T 



" dr 



dr 



(A3) 



L dr^ dr^ po \ dr / po \ dr / 



(A4) 



In the above two equations, ei„, eip, e2n and e2p are the functionals of nucleon densities /o„(r) and Pp{r) whose specific 
expressions are given in Ref. | [lO|] . Using the approximation of /i(a;) k, 2/x and l2{x) ~ 2 which are explained and 
employed in Ref. llQ], the following result can be obtained: 



In the simplified Myers-Swiatecki interaction, we have pU 



(A5) 



eir; 



Pn , Pp 

ai h au — 

Po Po 



f-ip — 



Pp , Pn 

ai h au — 

Po Po 



where a/_„ = ^(1 ± S,)a, thus the functional £q£ can be reduced to 



For the symmetric case with Pn ^ Pp ^ p/2, we have finally 



£2n = e2p = 0, 



£MS ^ ^|^^(v^)2 ^ _ «^)[(Vp„)2 + (Vpp)2]}. 

Po 



(A6) 



(A7) 



4if - ^aiVpr 
^Po 



(A8) 
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APPENDIX B: 



The formulas to calculate the interaction parameters from the nuclear matter quantities ai, Kq, J, L, and Kg will 
be given here for Skyrme, Myers-Swiatecki, and Tondeur interactions respectively. In Skyrme interaction, si, S2, ts, 
and a;3 can be calculated by the following equations: 



5l + i.2 



2 \2/3 5 

3tt^) 6(7-5) 



3(7-2) , 8701 -i^o 



T 



T 
Po 



(Bl) 



Sl + 2S2 = 



2 \2/3 



2(7 - 5) 



7 - 2 7(3 J -L) + Kg 



T 

5/3 ' 
Po 



(B2) 



i3 = 



16 



(7-5)(7-3) 



- 15ai 9 
f 5 



T 



(B3) 



_ 3r- 5(3J- L) - A'^ 1 
~ 2 isTo - 15ai - |T 2' 

Having si, S2, and ^3, to can be calculated by (|30|). Finally, a;o can be calculated by 



l^ 2 



7-2- 



57 J - (7 + 2)L + Kg 



T 
Po' 



Myers-Swiatecki interaction parameters can be calculated as: 

Ko - lOai 

a = — , 

T 



B = 



7 = 1 



4B(l + 7) 



5 if — 6ai 
18 T ' 

5 Kq — 15ai 
~ 9 T ' 



5B-J J 2B-J L 

+ 



a{4B + j)l B{l+j)T 25(1 + 7) T 



(B4) 

(B5) 

(B6) 
(B7) 
(B8) 
(B9) 



<^=i_ 2(1 + 7) 



3 3(45 + 7) 



1 - 



3 3J-L 



1 + 7 T 



Tondeur interaction parameters are 



37-2/^ 5 7 ai\T 
57-3V 37-2r/po' 



3 1 



(BIO) 

(Bll) 

(B12) 



2/3 ■ 

Po 



(B13) 
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TABLE I. The coefficients of volume energy ai, symmetry energy J, incompressibility Ko, density symmetry L and symmetry 
incompressibility Ks calculated from the various Skyrme interactions [Q(lst to 5th row), Myers-Swiatecki interaction [^(6th 
row), and Tondeur interaction [^(7th row), by original parameters, all in MeV. 7 is a model parameter in Eqs. § and (0) 
for Skyrme and Tondeur interactions respectively. EC is the equilibrium criterion calculated from the lefthand side of Eq. 
(^l|). As a comparison, the last two rows (labeled by CWS) present the result obtained by fitting these quantities directly to 
nuclear masses [Hj. ro is the nuclear radius constant in fm. 



EOS 


EC 


ro 




7 


ai 


Ko 


J 


L 


Ks 


SIII 


0.00080 


1.180 




6 


15.86 


355.5 


28.16 


9.88 


-393.9 


Ska 


-0.00001 


1.154 




4 


15.99 


263.1 


32.91 


74.62 


-78.45 


SkM 


0.00004 


1.142 




7/2 


15.77 


216.6 


30.75 


49.34 


-148.8 


SkM* 


0.00004 


1.142 




7/2 


15.77 


216.6 


30.03 


45.78 


-155.9 


RATP 


0.00049 


1.143 




18/5 


16.05 


239.6 


29.26 


32.39 


-191.3 


M-S 


0.00001 


1.140 






16.24 


234.4 


32.65 


49.88 


-147.1 


Tondeur 


0.00043 


1.145 




4 


15.98 


235.8 


19.89 


39.78 


-39.78 


CWS 


0.00000 


1.140 




4 


15.98 


217.5 


28.50 


64.32 


-101.3 


CWS 


0.00000 


1.140 




5 


16.10 


237.9 


28.50 


63.93 


-114.2 




TABLE II. Input values used to readjust the interaction parameters, ro 


in fm, others in MeV. 




Force 


ro 




ai 




Ko 


J 




L 


Ks 


Skyrme 


1.140 




15.97 




236.07 


29.25 




58.50 


-67.92 


M-S 


1.140 




15.97 




236.07 


29.25 




58.50 


-67.92 


Tondeur 


1.140 




15.97 




236.07 


29.25 




58.50 


-58.50 
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TABLE III. Readjusted Skyrme interaction parameters to, t^, xq, Si, and S2. Input values are ro = 1.140/m, ai = 15.97MeV , 
Ko = 236.07Mey", J = 29.25MeV, L = 58.50Mel/, and Ks = -67.92MeV. 



Force SIII Ska SkM SkM* RATP 

7 6 4 7/2 7/2 18/5 

to{MeVfm^) -1405.521 -1792.320 -2372.518 -2372.518 -2179.119 

t3{MeVfrn<) -14402.55 12794.56 12584.33 12584.33 11940.95 

xo 0.06956 0.13735 0.19759 0.19759 0.18018 

X3 0.38368 0.38368 0.38368 0.38368 0.38368 

siiMeVfnC') 642.825 -42.389 71.813 71.813 55.499 

S2{MeVfm^) -473.186 84.802 -8.196 -8.196 5.090 



TABLE IV. Critical point [pc, 5c) predicted by Skyrme, Myers-Swiatecki, and Tondeur interactions respectively, ro in fm, pc 
in fm~^, and Cc in MeV. For each item, the first line is given by original interaction parameters, the second line by readjusted 
parameters shown in Table 3 for Skyrme interactions while by Eqs. ( |60| ) and (^l|) for Myers-Swiatecki and Tondeur interactions 
respectively. 





SIII 


Ska 


SkM 


SkM* 


RATP 


M-S 


Tondeur 


ro 


1.180 


1.154 


1.142 


1.142 


1.143 


1.140 


1.145 




1.140 


1.140 


1.140 


1.140 


1.140 


1.140 


1.140 


Sc 


0.8385 


0.8647 


0.8390 


0.8421 


0.8303 


0.8213 


0.8732 




0.8772 


0.8980 


0.8908 


0.8908 


0.8920 


0.8988 


0.7697 


Pc 


0.07173 


0.02416 


0.02345 


0.02420 


0.03892 


0.03039 


0.03081 




0.02732 


0.02969 


0.02825 


0.02825 


0.02851 


0.02643 


0.03131 


Cc 


3.9019 


1.5852 


1.2572 


1.2814 


1.9898 


1.1031 


2.6142 




1.8894 


1.8505 


1.7025 


1.7025 


1.7311 


1.1280 


2.6304 
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FIG. 1. (a) Normalized energy per nucleon e/ai versus relative nucleon density p/po calculated by original interaction 
parameters at 5 = 0, (b) e versus p/po calculated by readjusted interaction parameters at J = of the various Skyrme 
interactions (solid lines), Myers-Swiatecki interaction (dot-dashed line), and Tondeur interaction (dashed line). The solid lines 
in the righthand side of the plot in (a) from top to bottom correspond to SIII, Ska, RATP, SkM, and SkM* interactions 
respectively. The difference between SkM and SkM* is negligible and Myers-Swiatecki is almost coincident with RATP. 

FIG. 2. (a) Normalized nuclear incompressibility K/ Kq versus relative nucleon density p/po calculated by original interaction 
parameters at 5 = 0, (b) K versus p/po calculated by readjusted interaction parameters at 5 = of the various Skyrme interac- 
tions (solid lines), Myers-Swiatecki interaction (dot-dashed line), and Tondeur interaction (dashed line). On the righthand side 
of the plot in (a) the solid lines from top to bottom correspond to SIII, Ska, RATP, SkM, and SkM* interactions respectively. 
The difference between SkM and SkM* is negligible. In (b), the solid lines from top to bottom are due to SkM, SkM*, RATP, 
Ska, and SIII, where SkM* is identical to SkM. Tondeur is coincident with Ska; SkM and RATP almost overlap. 

FIG. 3. (a) Normalized symmetry energy J(p)/J versus the relative nucleon density p/po calculated by original interaction 
parameters, (b) J(p) versus p/po calculated by readjusted interaction parameters of the various Skyrme interactions (solid 
lines), Myers-Swiatecki interaction (dot-dashed line), and Tondeur interaction (dashed line). On the righthand side of the plot, 
in (a), the solid lines, from top to bottom, correspond to Ska, SkM, SkM*, RATP, and SIII interactions respectively. The 
difference between SkM and SkM* for p/po < 1.3 is negligible. In addition, the Myers-swiatecki's is almost coincident with 
that of SkM's. In (b), the difference among these curves is negligible. 

FIG. 4. Equilibrium density pm as function of the relative neutron excess 5, calculated by (a) original parameters and 
(b) readjusted parameters of the various Skyrme interactions (solid lines), Myers-Swiatecki interaction (dot-dashed line), and 
Tondeur interaction (dashed line). The solid lines, in the middle range of 5 from top to bottom in (a) correspond to SIII, 
RATP, SkM, SkM*, and Ska interactions respectively, the difference between SkM and SkM* is very small. In (b) the solid 
lines from top to bottom correspond to Ska, RATP, SkM, SkM*, and SIII interactions respectively. SkM and SkM* are the 
same whereas RATP, SkM and Myers-Swiatecki almost overlap. 

FIG. 5. Equilibrium energy per nucleon em as function of the relative neutron excess 5, calculated by (a) original parameters 
and (b) readjusted parameters of the various Skyrme interactions (solid lines), Myers-Swiatecki interaction (dot-dashed line), 
and Tondeur interaction (dashed line). The solid lines, in the righthand side, in (a) correspond to RATP, SkM, SkM*, SIII, and 
Ska interactions from top to bottom respectively. SIII, Ska and Tondeur almost overlap, whereas the difference between SkM 
and SkM* is negligible. In (b) the solid lines from top to bottom correspond to SIII, SkM, SkM*, RATP, and Ska interactions, 
respectively, where SkM and SkM* are the same. SkM, RATP and Ska are almost coincident. 

FIG. 6. Equilibrium incompressibility Km as function of the relative neutron excess 5, calculated by (a) original parameters 
and (b) readjusted parameters of the various Skyrme interactions (solid lines), Myers-Swiatecki interaction (dot-dashed line), 
and Tondeur interaction (dashed line). The solid lines in (a) correspond to SIII, RATP, Ska, SkM*, and SkM interactions from 
top to bottom in the middle range of 5, respectively. Ska, SkM* and SkM almost overlap. In (b) the solid lines from top to 
bottom correspond to Ska, RATP, SkM, SkM*, and SIII interactions respectively; SkM and SkM* are the same. 

FIG. 7. Neutron matter EOS e(p, 1) versus nucleon density p, calculated by (a) original parameters and (b) readjusted pa- 
rameters of the various Skyrme interactions (solid lines), Myers-Swiatecki interaction (dot-dashed line), and Tondeur interaction 
(dashed line). Black diamonds denote the data taken from Ref. The solid lines, in the right hand side of (a), correspond 
to Ska, SkM, SkM*, RATP, and SIII from top to bottom. In (b), the solid lines, from top to bottom, correspond to SIII, SkM, 
SkM*, RATP, and Ska, where SkM and SkM* are identical. 



FIG. 8. Symmetric nuclear matter EOS e(p, 0), calculated by (a) original parameters and (b) readjusted parameters of the 
various Skyrme interactions (solid lines), Myers-Swiatecki interaction (dot-dashed line), and Tondeur interaction (dashed line). 
The full dots stand for the results taken from Ref. the crosses denote the results taken from Ref. [^. The solid lines, in 
the right hand side of (a), correspond to SIII, Ska, RATP, SkM, and SkM*, from top to bottom, where SkM and SkM* are 
identical. In (b), the solid lines, from top to bottom, correspond to SkM, SkM*, RATP, Ska, and SIII, where SkM and SkM* 
are identical. The Tondeur's is very close to the Ska's. 
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FIG. 9. Neutron matter EOS e(p, 1), calculated by (a) original parameters and (b) readjusted parameters of the various 
Skyrme interactions (solid lines), Myers-Swiatecki interaction (dot-dashed line), and Tondeur interaction (dashed line). The 
full dots stand for the results taken from Ref. Q, the crosses denote the results taken from Ref. j^. The solid lines, in the 
right hand side of (a), correspond to Ska, RATP, SkM, SkM*, and SIII, from top to bottom. In (b), the solid lines, from top 
to bottom, correspond to SIII, SkM, SkM*, RATP, and Ska, where SkM and SkM* are identical. 
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